An Inflationary Non-singular Quantum Cosmological Model 
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A stiff matter-dominated universe modeled by a free massless scalar field minimally coupled 
to gravity in a Friedmann-Lemaitre-Robertson- Walker (FLRW) geometry is quantized. General- 
ized complex-width gaussian superpositions of the solutions of the Wheeler- De Witt equation are 
constructed and the Bohm-de Broglie interpretation of quantum cosmology is applied. A planar 
dynamical system is found in which a diversity of quantum bohmian trajectories are obtained and 
discussed. One class of solutions represents non-singular inflationary models starting at infinity past 
from fiat space-time with Planckian size spacelike hypersurfaces, which inflates without inflaton but 
due to a quantum cosmological effect, until it makes an analytical graceful exit from this inflationary 
epoch to a decelerated classical stiff matter expansion phase. 
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I. INTRODUCTION 

For more than 25 years, inflation [l[ has been consid- 
ered a paradigm to solve, at the same time, standard 
cosmological puzzles related to initial conditions like the 
flatness, horizon, and isotropy problems, and, as a bonus, 
astroparticle issues like the monopole excess. More im- 
portant, it also predicts that primordial fluctuations, as- 
sumed to be of quantum origin, could be enhanced to 
the level required to trigger large scale structure forma- 
tion, with an almost scale-invariant spectrum [2J, which 
is confirmed by observations 0|. 

The inflation paradigm is also endowed with two spe- 
cific problems, that may ultimately be related, namely 
the meaning of the trans-Planckian [4| perturbations and 
the existence of a past singularity [5j]. Concerning the 
latter, the existence of an initial singularity is one of the 
major drawbacks of classical cosmology. In spite of the 
fact that the standard cosmological model, based in the 
classical general relativity theory, has been successfully 
tested until the nucleosynthesis era (around t ~ Is), the 
extrapolation of this model to higher energies leads to a 
breakdown of the geometry in a finite cosmic time. This 
breakdown of the geometry may indicate that the clas- 
sical theory must be replaced by a quantum theory of 
gravitation: quantum effects may avoid the presence of 
the singularity, leading to a complete regular cosmologi- 
cal model. 
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Among the fundamental questions that come from the 
quantization of the universe as a whole, one of the most 
important concerns the interpretation of the wave func- 
tion of the Universe. In order to extract predictions from 
it, the Bohm-de BrogliefBdB) ontological interpretation 
of quantum mechanics 0, 0] has been proposed [8|, [^, [Io| , 
since it avoids many conceptual difficulties that follow 
from the application of the standard Copenhagen inter- 
pretation to a unique system that contains everything. 
In opposition to the latter one, the ontological interpre- 
tation does not need a classical domain outside the quan- 
tized system to generate the physical facts out of poten- 
tialities (the facts are there ab initio because the posi- 
tions and trajectories of the particles (called bohmian 
trajectories) are considered to be part of objective re- 
ality, and hence it can be applied to the Universe as a 
whole. There are other alternative interpretations which 
can be used in quantum cosmology, like the many worlds 
interpretation of quantum mechanics [TlT ] , but we will not 
consider them here because they are probabilistic inter- 
pretations in essence. As we know [12j . it is very dif- 
ficult to obtain from the the Wheeler-DeWitt equation, 
when applied to a closed universe, a probabilistic inter- 
pretation for their solutions because of its hyperbolic na- 
ture (see however further approachs [HI)- In the case 
of the Bohm-de Broglie interpretation, probabilities are 
useful but not essential, as long as objective trajectories 
(universe histories) can be calculated and their proper- 
ties studied. Probabilities can be recovered, as it has 
been suggested many times [HI, EH EH, E3 , at the semi- 
classical level, where a probability measure can be con- 
structed with the quantum solutions. Hence, we can take 
the Wheeler-DeWitt equation as it is, without imposing 
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any probabilistic interpretation at the most fundamental 
level, but still obtaining information using the Bohm- 
de Broglie interpretation, and then recover probabilities 
when we reach the semiclassical level. With this inter- 
pretation in hands, one can ask if the quantum scenario 
predicted by the wave function of the Universe is free of 
singularities, and which type of classical universe emerges 
from the quantum phase. 

Quantum cosmology, in this framework, exhibits 
bouncing solutions [18[ which can be interpreted as truly 
avoiding the singularity. Some of these bouncing models 
provide solutions to many cosmological puzzles, and may 
also yield a scale invariant cosmological perturbations as 
in inflationary models [2(j. Bouncing models are also 
obtained in loop quantum cosmology [l9(. 

Some types of bouncing models were obtained in 
[2H . |22| . In these references, the matter content of the 
primordial model was considered to be stiff matter , 
modeled by a free massless scalar field minimally coupled 
to gravity. The Wheeler-DeWitt equation turns out to 
be a two-dimensional Klein-Gordon equation, and gaus- 
sian superpositions, with positive width, of negative and 
positive frequency modes solutions were considered. In 
the flat case, a two dimensional dynamical system for 
the bohmian trajectories was obtained, yielding a vari- 
ety of possibilities: big bang-big crunch models, oscil- 
lating universes, bouncing solutions, and ever expanding 
(contracting) big bang (big crunch) models with a period 
of acceleration in the middle of their evolution. No non- 
singular purely expanding model with a primordial ac- 
celerated phase (a non-singular inflationary model) was 
obtained (models of this type were obtained only with 
a non-minimal coupling between the scalar and gravita- 
tional fields, [H). 

In the present paper we g eneralize the gaussians su- 
perpositions of Refs. [2l|, |22( | to gaussians with complex 
widths and non negative real part. We obtain a richer 
two dimensional dynamical system. When the real part 
of the width is made zero, we obtain non-singular infla- 
tionary models which expand accelerately in the infinity 
past from flat universes with finite (not zero) size spatial 
sections, and is smoothly connected to a classical decel- 
erated stiff matter expanding phase. It has features of 
the pre-big bang model [13] and the emergent model [28| 
without graccfull exit problem. 

The paper is organized as follows: in the next section 
the classical model is presented. Section III is devoted to 
its quantization and the corresponding Wheeler-DeWitt 
equation is obtained. Generalized gaussian superposi- 
tions of the quantum solutions, and their corresponding 
dynamical system are studied in section IV. In Section 
V, the non-singular inflationary model is studied and its 
properties are presented and discussed. In section VI we 



present our conclusions. 



II. THE CLASSICAL MINISUPERSPACE 
MODEL 

The model we take contains a massless free scalar field 
(stiff matter), and the total lagrangian reads 



(1) 



where we are using natural units h = c = 1, and I 2 = 
87rG/3, which is the Planck length squared in these units. 

We will consider the spatially homogeneous and 
isotropic space-time line element, 



ds 2 = -N 2 dt 2 



a(tY 



(l + er 2 /4) ; 



: [dr 2 +r 2 (d9 2 +sm 2 (d)dcp 2 )} 



(2) 

where the spatial curvature e takes the values 0, 1,-1. 
Inserting this line element into the lagrangian ([T]), and 
omitting a total time derivative, we obtain the following 
minisuperspace action: 



S = 



-a z aV NeaV (j> 2 a 3 V 



Nl 2 



P 



2N 



dt 



(3) 



where V is the total volume divided by a 3 of the space- 
like hypersurfaces, which are supposed to be closed. V 
depends on the value of e and on the topology of the hy- 
persurfaces. For e = 0, V can have any value because the 
fundamental polyhedra of e — hypersurfaces can have 
arbitrary size (sec Rcf. 29]). In the case of e = 1 and 
topology S 3 , V = 2tt 2 . 

Usually, the scale factor has dimensions of length be- 
cause we use angular coordinates in closed spaces. Also, 
in natural units, the scalar field has dimensions of the 
inverse of a length. Hence we will define the dimensio- 
less quantities a = \/2Va/l, <j> = l<p/\/2. Calculating the 
hamiltonian, and omitting the bars, yields, 



H = 



I 



2a 2a 3 



at 



(4) 



As y2V jl appears as an overall multiplicative constant 
in the hamiltonian, we can set it equal to one without 
any loss of generality, keeping in mind that the physical 
scale factor which appears in the metric is la/y/2, not 
a. We can further simplify the hamiltonian by defining 
a = ln(a) obtaining 



H 



N 



2 exp(3a) 



- Pi + Pl - e exp(4a) 



, (5) 



where 



1 A perfect fluid with a "super-ri gid" equation of state p - 
was proposed long time ago in 12511 . 



: p, that 



Pa 



e 3a a 
N 



e ia 4> 
N 



(6) 
(7) 
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The momentum p^ is a constant of motion which we will 
call k. 

The classical solutions are, in the gauge N = 1 (cosmic 
time): 



1) For e = 0: 



±ci + ci 



(8) 



where c\ is an integration constant. In terms of cosmic 
time r they read: 



a = e a = 3kr 1/3 , 



C2 



(9) 
(10) 



III. QUANTIZATION AND THE BOHM-DE 
BROGLIE INTERPRETATION 



Let us now quantize the model. The Wheeler-DeWitt 
equation is obtained through the Dirac quantization pro- 
cedure, where the wave function must be annihilated by 
the operator version of the Hamiltonian constraint. For 
the case of homogeneous minisuperspace models, which 
have a finite number of degrees of freedom, the minisu- 
perspace Wheeler-De Witt equation reads 



W(p",« M )*(g) = 



(15) 



The quantities p^,^ are the phase space operators re- 
lated to the homogeneous degrees of freedom of the 
model. Usually this equation can be written as 



The solutions contract or expand forever from a singular- 
ity, depending on the sign of k, without any inflationary 
epoch. 

2) For e = 1: 



k 



cosh(20 — ci) ' 



(11) 



where c\ is an integration constant, and from the conser- 
vation of Ptf, we get 



1 (q) dm 

2 dq p dq<j 



+ UfoMq) = 



(16) 



where f pa (qp) is the minisuperspace DeWitt metric of 
the model, whose inverse is denoted by f pa {q P )- 

Writing iff in polar form, iff = Rexp(iS), and substi- 
tuting it into (fTr?|) . we obtain the following equations: 



(17) 



(12) 



The cosmic time dependence is complicated and we will 
not write it here. These solutions describe universes ex- 
panding from a singularity till a maximum size and con- 
tracting again to a big crunch. Near the singularity, these 
solutions behave as in the flat case. There is no inflation. 



3) For e = 1: 



k 



sinh(2(/> — c\) 



(13) 



where c\ is an integration constant, and again, from the 
conservation of p$ we get 



k = e ia (j) 



(14) 



As before, the cosmic time dependence is complicated 
and we will not write it here. These solutions describe 
universes contracting forever to or expanding forever 
from a singularity. Near the singularity, these solutions 
behave as in the flat case. There is no inflation. 

Hence, in all models there is at least one singularity 
and no acceleration phase, as it should be for a classical 
stiff matter fluid. 



fpcilp) 



'dq p 



R" 



dS 
dq a 



= 



(18) 



where 



Q(qu) 



i 



d 2 R 



2R Jpu dq p dq a 



(19) 



is called the quantum potential. 

The Bohm -de Broglie interpretation applied to Quan- 
tum Cosmology states that the trajectories q^t) are real 
independently of any observations. Equation (JT7J) repre- 
sents their Hamilton-Jacobi equation, which is the clas- 
sical one added with a quantum potential term Eq. (|19[ 
responsible for the quantum effects. This suggests to de- 
fine 



dS 



(20) 



where the momenta are related to the velocities in the 
usual way: 



P p = f 



pa 



1 dqcr 
N dt ' 



(21) 
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To obtain the quantum trajectories we have to solve 
the following system of first order differential equations, 
called the guidance relations: 



9S(q, 
dq P 



- = f Pa T-r<l°- 



1 

N~- 



(22) 



Eqs. (|22j) arc invariant under time reparamctrization. 
Hence, even at the quantum level, different choices of 
N(t) yield the same space-time geometry for a given non- 
classical solution q a (t). 

There is no problem of time in the Bohm-de Broglie in- 
terpretation of minisuperspace quantum cosmology [30| . 
This is not the case, however, for the full supcrspace, 
see [3l| , although the theory remains consistent, see 

Mm)- 

Let us then apply this interpretation to our minisuper- 
space model. The operator version of Eq. (O, with the 
factor ordering which makes it covariant through field 
redefinitions, reads 



1 



<9 2 * a 2 * 



2e 3 » V 9a 2 d<j) 2 



ee 4a * 



(23) 



from where one can read that f pa = r/ pa e , and r\ pa 
is the Minkowski metric in two dimensions. Comparing 
Eq. (J23J) with Eqs. (31} and ([T7J) yields for the quantum 
potential 



Q(a, ( 



„3a 



e"" d 2 R 
R pa dq p dq a 



I 

R 



d 2 R d 2 R 
da 2 dd> 2 



The guidance relations ([22]) read 



dS_ 

da 



e 3a a 
N 



(24) 



(25) 



yielding, 



dudv 



(28) 



The solutions are: 
1) For e = 0: 

In this case the general solution is 



V(u,v) = F(u) + G(v) , (29) 
where F and G are arbitrary functions. Using a separa- 
tion of variable method, one can write these solutions as 
Fourier transforms given by 



%(u,v) = J dkll(k)exp(iku) + J dkV (k) exp(ikv) 

U and V also being arbitrary. 

Writing the solution (|2"9"|) in polar form, 



(30) 



= R+e lb + +R-e 



JS- 



where 



R+ = R{u) S+ = S{u) 
R- = R(v) S- = S(v) 



dS_ e^6 
90 ~~ N 

We can write this equation in null coordinates, 



u = —=(a + 

v = —= (a - 
V2 



1 / 

a = — {u + v) 



(26) 



one obtains 



R = \IR 2 + + R 2 _ +2R+R-cos(S+ - S-) 
R+ sin(S + ) + R-sin(S-) 



S = arctan 



R+ cos(S + ) + R- cos(S_) 



(27) The derivative of S with respect to some variable x reads 



OS = + R^ + + ^r) R + R- cos (S + s-) + (R- d -§±- R + 9 4r) sin (S + flL) 

dx 



R 2 + + R 2 _+ 2R+R_ cos{S+ - S-) 



These equations will be used in the next section. 
2) For e ^ 0: 



In this case the general solution reads 



v) 



dk U (k) exp 



2V2 



4V2fc 



dk V(k) exp 



2\/2 



e 2 ^") , (31) 



4\/2fc 
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where / and g are arbitrary functions. In Ref.[2l| these 
solutions were expanded in terms of Besscl functions. 



IV. GENERALIZED GAUSSIAN 
SUPERPOSITIONS 

In Ref. [U, [22| we made gaussian superpositions of the 
solutions with the choice U(k) = V(±k) — A(k), with 
A{k) given by 



A(k) = cxp 



(k - V2d) 2 



(32) 



with a 2 > 0, and the presence of y/2 above is just for fur- 
ther convenience. Bouncing non-singular solutions were 



obtained in [2l|, and expanding singular models with a 
period of acceleration between decelerated phases were 
proposed in [22I ]. 

In this paper we will consider the more general case 
where the parameter a 2 in (|32[) is given by a complex 
number: a 2 + i4h, where h is an arbitrary real number. 
Under this asumption we have: 



A(k) = cxp 



(k - V2d) 2 
a 2 + i4h 



(33) 



From now on we will consider only flat spatial sections. 

Integrating with U(k) = V(k) = A(k), we obtain 
the solution: 



&(u,v) = y/a 4 + 16h 2 



i arctan . / + — ct_ 

e V V" + 16h +" ^ exp 



—u 2 + i(-hu 2 + V2di 
4 



exp 



-v 2 + i(-hv 2 + V2dv) 



(34) 



To obtain the quantum trajectories, we have to calculate the phase S of the above wave function and substitute it 
into the guidance equations. Computing the phase and recovering the original variables a, (ft, we have 



_ / / \/c7 4 + 16h 2 — <j 2 \ h o 2 f / CT 2 Q(/^ 

b = arctan \ / — I + da (a + (ft ) + arctan< tanh ] ta,n[(ft(ha — d)\ 

Vy Vc 4 + I6h 2 + a 2 ) 2 \ 4 

which, after substitution in Eqs (I25l26p . yields a planar system given by: 



(35) 



and 



j\t r a 2 <fi sm[2(j)(ha — d)] + 4h(f) smh . 

4(d - ha) + 



4e 3a 



N 

4 e 3a 



Ah(ft + 



cosh (^2) + cos[24>{ha - d)} 
a 2 a sm[2<f)(ha - d)} + A(ha - d) sinh 



cosh 



(^) + cos[2(ft(ha - d)} 



The norm of the solution (|34[) is given by 



=: f{ot,4>), 



g{a,4>). 



(36) 



(37) 



R = V2^ </a A + 16h 2 e -V(« +4> ) J cosh 



cos[2(j)(ha — d)] 



(38) 



Equations (|3l))) . (|3T|) give the directions of the geo- 
metrical tangents to the trajectories which solves this 
planar system. We shall work in the gauge N = 
e 3a . By plotting the tangent direction field, it is pos- 
sible to obtain the trajectories. Due to the symme- 



tries f(a,(ft;h,d) = /(— a, — </>; — h, d), g(a, — (ft; h, d) = 
g(-a,-(ft;-h,d), and f(a,(ft;h,d) = -f(-a,<ft;h,-d), 
g(a,(ft;h,d) = g{— a, (ft; h, — d), one concludes that a 
change in sign of h corresponds to a inversion around 
the origin with time reversion, and a change in sign of 
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FIG. 1: Field plot for the system d36j , dHZJ) , giving the di- 
rection of the geometrical tangent to the trajectories, for the 
values a 2 = 2, h = 1/8, and d = —1. 
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FIG. 2: Field plot for the system (|36} , (|37|) , giving the di- 
rection of the geometrical tangent to the trajectories, for 
a 2 — 2, h = 0.5, and d = — 1. For a bigger value of h the cen- 



ter points are farther from the axis a 
in the sign of a when <j> is big enough. 



0. Note the change 



d corresponds to a reflexion in the (f> axis. Hence, one 
can make definite choices of sign for these parameters 
without loss of generality. 

Field plots of this planar system arc shown in 
Figs (|H2l3l4p . for the choice of parameters {a 2 = 2, h = 
1/8}, {a 2 = 2,h = 0.5}, and {a 2 = 2, h = 5} (two por- 
traits), respectively, all with d = — 1. 

The line 4> = divides configuration space in two sym- 
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FIG. 3: Field plot for the planar system (J36J) (J3TJ) for cr 2 = 
2,h — 5, and d = — 1. Two trajectories are depicted: one 
representing a bouncing universe spending a long time on the 
bounce and other which corresponds to a universe which be- 
gins and ends in singular states ( "big bang - big crunch" uni- 
verse). Note the change in the sign of q when cj> is big enough. 



metric regions as long as /(a, <j)) = /(a, — </>) in Eq. (|36p . 
and g(a,(f>) = —g(a,—<p) in Eq. (|37ll . This can be seen 
in Fig J!]). 

The line a = contains all the nodes of this system, as 
it is shown in all figures. They appear when the denom- 
inator of the above equations, which is proportional to 
the norm of the wave function (see Eq. (|38p). is zero. No 
trajectory can pass through these points. They happen 
when a = and cos(2d(f>) = —1, or cf) = (2n + l)ir/(2d), 
n an integer, with separation ir/d. 

The center points appear when the numerators are 
zero, their locations depend on the values of h, a 2 and 
d, and they are not on the line a = 0, unless h = (case 
of Rcfs. [2l|, HI]). Note in Fig {I), where h is relatively 
much smaller than a 2 , that the centers are close but not 
on the line a — 0. 

Apart the changes in the sign of d, which happen 
around the center points, there are regions with different 
signs of a when \<f>\ » 0, \a\ » where the hyperbolic 
functions dominate over the trigonometric. For instance, 
considering (f> > 0, if d < 0,h > 0, and a > for a fixed, 
Eq. ([36]) shows that a > for small (j> when the cosh term 
dominates, but it change signs when cf> is big enough in 
order for the sinh term be greater than the cosh term in 
this equation, changing the sign of a. This situation is 
depicted in Figs (|2|3|4[) . 

Note finally that the classical solutions for a(r) 
(a(r) cx t 1 / 3 ) and 0(r) are recovered when | a |— > oo 
or | (j) | — ► oo , and none of them are null. 

We can see plenty of different trajectories, depending 
on the initial conditions and parameter values, in the 



FIG. 4: Field plot for the planar system (|36|) . (|37|) giving the 
direction of the geometrical tangent to the trajectories for 
a 2 = 2, h — 5, and d = — 1. The trajectory for a universe 
coming from a singularity, experiencing a long static phase 
and finally expanding, is depicted. 



Figs (|H2l3l4p . Near the center points we can have os- 
cillating universes without singularities. One may have 
universes expanding classically from a singularity, experi- 
encing quantum effects in the middle of their expansion, 
with possible accelerating and/or static phases, and re- 
covering their classical expansion behaviour for large val- 
ues of a. 

There are big bang-big crunch models, and bouncing 
models, where the bounce may take long and connects 
two asymptotically classical contracting and expanding 
phases (see Fig 

Note that, if we choose V(k) = A(— k) (see Eqs 
(|30I32[) ). one obtains the same field plots as above with 
the axis a and <fi interchanged, obtaining more possibili- 
ties for bouncing models. 

Finally, there is the special situation of a 2 = 0, which 
leads to qualitative different solutions, as we will see in 
the next section. 



V. NON-SINGULAR INFLATIONARY 
BOHMIAN TRAJECTORIES 

There is an interesting case which is obtained when 



(j = in the gaussian ([33 
Hence, A(k) reads 

A(k) = exp 



(k - V2d) 2 
1 Ah 



(39) 



Then the wave function (|3"4"|) reduces to 



FIG. 5: Field plot for the linear planar system (|42|) (|43j) for 
a 2 = 0, h — 0.5, and d = — 1. A class of trajectories represent 
an universe that begins to inflate quantum mechanically from 
a Planckian size closed flat space-time in the infinite past, 
and when it becomes large enough it undergoes an analytical 
graceful exit to a decelerated classical stiff matter expansion 
phase. 



V(u,v) = 2y/w\h\ 



exp i (^-hu 2 + \/2du + — ^ 



+ exp i ( — hv 2 + \/2dv + ^ 



(40) 



Its norm is given by R = 4y7r[/ij cos[4>(ha — d)], yield- 
ing the quantum potential 



(ha — d) 2 - h 



2 J,2 



(41) 



The guidance relations given by (|36|) and (|37|) now re- 
duce to 



ha — d 



4> = —hcj> 



(42) 



(43) 



The equations (f4"2"| and (|13")) represent a linear dynamical 
system. The only critical point (<fr = 0, a = 4) is a saddle 
point and, as it is well known, it represents an unstable 
equilibrium. 

The field plot of these solutions is depicted in Fig ([5]) 
for d = —l,h = 0.5. Note that there are two regions of 
different signs of a separated by the line a = d/h. 

In this case there are analytical solutions, which read 



e d l h exp(a e' lt ) and 4> = </> e 



-ht 



(44) 



where ao and 4>o are integration constants, and remem- 
bering that the time parameter t is related to cosmic time 
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r through t = J dte 3a(t) ^t-t = Ei(3a e ht )/h, where 
Ei(x) is the exponential-integral function. 

These solutions represent ever expanding or contract- 
ing non-singular models, depending on the sign of h. Let 
us consider the physical situation of expanding solutions 
with h > 0. The Hubble and deceleration parameters 
a' / a and a" /a read (a prime denotes a derivative in cos- 
mic time) 



and 



a 



a 
a 

a h' 



a he 



III 



e ht {l-2a e ht ) 



and the scalar curvature is 



n = - 



6a /i 2 M,-, ht\ 

n u 



(45) 



(46) 



(47) 



There are three important phases in this model. For 
t << the Universe expands accelerately from its min- 
imum size ciq = e d / h (remember that for the physical 
scale factor one has ciQ hys = le d l h /\J2V) 7 which occurs 
in the infinity past t — ► — oo when the curvature is 
null but increasing while scale factor grows. The scalar 
field is very large in that phase. For t >> the Uni- 
verse expands decelerately, the scale factor is immensely 
big, the scalar field becomes negligible and the curva- 
ture approaches zero again. The transition occurs when 

Wtran = -ln(2a ). 

Around ht = one has 

a a e aa+d ' h [l + a ht + (a Q h 2 + a 2 a h 2 )t 2 /2! + ...]. (48) 

If ao >> 1 (and hence t > ft ran , which means in the de- 
celeration phase), one can write a ~ e ao+d / h cxp(aoht). 
In that case, from r = J dta 3 (t), one obtains that 
a oc (r — To) 1 / 3 and <fi' cx 1/r cx 1/a 3 , as in the clas- 
sical regime. 

Collecting all these phases, and considering a >> 1, 
one has a non-singular ever expanding Universe, starting 
with a constant and finite size (which may be of the or- 
der of the Planck length if d = 0, a Planckian flat space) 
in the infinity past, which inflates afterwards till it at- 
tains an almost 2 classical decelerating expanding regime, 
with a size e a ° times bigger than it was initially. At this 
time, radiation may start to dominate over the scalar 
field since it becomes much smaller than in the infinity 
past, and its energy density goes may like 1/a 6 . There 
are no event nor particle horizons. Hence we have an 



2 We would like to say that it is not necessary that classical be- 
haviour appears when the Universe is large, see |23H . In some 
circumstances, a quantum behavior should be desirable, as in 
[22l | to produce late acceleration in the Universe, or to avoid a 
Big Rip, see [3 . 



inflationary non-singular model which can be smoothly 
joined to the standard model. 

Some types of bouncing model were obtained in [2lj . 
In this reference, the matter content of the primordial 
model was considered to be stiff matter modeled by a 
free massless scalar field minimally coupled to gravity. 
The Wheeler-DeWitt equation turns out to be a two- 
dimensional Klein-Gordon equation, and gaussian su- 
perpositions, with positive width, of negative and pos- 
itive frequency modes solutions were considered. In the 
flat case, a two dimensional dynamical system for the 
bohmian trajectories was obtained, yielding a variety 
of possibilities: big bang-big crunch models, oscillating 
universes, bouncing solutions, and ever expanding (con- 
tracting) big bang (big crunch) models with a period of 
acceleration in the middle of their evolutions. No non- 
singular purely expanding model with a primordial ac- 
celerated phase (a non-singular inflationary model) was 
obtained (models of this type were obtained only with 
a non-minimal coupling between the scalar and gravita- 
tional fields, (Hi). 



VI. CONCLUSION 

In the present paper we generalized the gaussians su- 
perpositions of Refs. [2l| [22j to gaussians with complex 
widths with non negative real part. We obtained a richer 
two dimensional dynamical system, with oscillating uni- 
verses without singularities, universes expanding classi- 
cally from a singularity, experiencing quantum effects 
in the middle of their expansion with possible acceler- 
ating and/or static phases whose duration depend on 
the three free parameters of the gaussian (its complex 
width and the location of its center) and two initial con- 
ditions, which recover their classical expansion behaviour 
for large scale factors. There are also big bang-big crunch 
models, and bouncing models, where the bounce may 
take long, and connects two asymptotically classical con- 
tracting and expanding phases. 

The most interesting solutions occur when the real part 
of the width is made zero. Then it was obtained an- 
alytical non-singular inflationary models which expand 
accelerately in the infinity past from flat universes with 
finite (not zero) size spatial sections, and are smoothly 
connected to a classical decelerated stiff matter expand- 
ing phase. 



It is like the pre-big bang model [27J, with 
a minimum volume spatial section in the infinity past 
(which can be of the order of the Plank volume), or as 
the emergent model [28| for flat spatial sections, with- 
out any gracefull exit problem. Sufficient inflation can 
be obtained with reasonable choices of initial conditions 
(ao > 70). The picture is then of a Universe that be- 
gins to inflate quantum mechanically from a Planckian 
size closed flat space-time in the infinity past, and when 
it becomes large enough it makes a smooth transition to 
a decelerated classical stiff matter expansion, which will 
be eventually be dominated by radiation before nuclc- 
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osynthcsis. Some reheating process must then also be 
presented in our model, and this could happen through 
interactions of the type \J— gvip^^J, where \& represents 
some fermion field, v is a coupling constant, and the 
scalar field decays into these fermions. Such fermions 
would not be present in the beginning but they would 
have been produced when the universe got bigger and 
transitatcd to the decelerated phase. Such extra term 
would not be relevant in the hamiltonian which describes 
the model when the universe was small (it is proprtional 
to a 3 ) and empty of these fermions, justifying our quan- 
tum solutions in this regime. However, when the fermions 
begin to be produced afterwards, our quantum solutions 
could not be anymore reliable. At this stage, when the 
model becomes decelerated and gets closer to classical 
behaviour, this solution must be connected to the classi- 
cal radiation plus stiff matter solution . The nonclassical 
behaviour of our bohmian quantum solution for t > 
should then be of course discarded, they are out of the 
range of the validity of our assumptions. Therefore, our 
non-singular inflationary solution should be utilized only 
as a zeroth-order approximation to a realistic cosmolog- 
ical inflationary model, in the same way as the de Sitter 



model or power law solutions are used for usual inflation 
[32| . The improvement of the model, including reheating 
and a radiation field deserves further research as there 
are many possibilities to be investigated. Another im- 
portant step forward is to evaluate the evolution of cos- 
mological perturbations using the formalism developpcd 
in Ref. [2(|, and confront it with observations [||. This 
will be the subject of our future publications. 
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